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Abstract 

A critical concern in the fabrication of targets for inertial confinement fusion (ICF) is ensuring that 
the hydrogenic (D 2 or DT) fuel layer maintains spherical symmetry. Solid layered targets have 
structural integrity, but lack the needed surface smoothness. Liquid targets are inherently smooth, 
but suffer from gravitationally induced sagging. One method to reduce the effective gravitational 
field environment is freefall insertion into the target chamber. Another method to counterbalance 
the gravitational force is to use an applied magnetic field combined with a gradient field to induce a 
magnetic dipole force on the liquid fuel layer. Based on time dependent calculations of the 
dynamics of the liquid fuel layer in microgravity environments, we show that it may be possible to 
produce a liquid layered ICF target that satisfies both smoothness and symmetry requirements. 

Introduction 

In the Inertial Confinement Fusion (ICF) concept, 1 multiple high power laser or ion beams are 
focused on a small spherical target containing condensed hydrogen fuel. The intense pressure 
generated by the incident beams cause the target to implode to high density (~1000x liquid density) 
and temperature (>100 million K) with the subsequent release of thermonuclear energy. For ICF 
to be economically attractive, the gain (thermonuclear energy out divided by beam driver energy 
in) must be large enough to overcome driver inefficiencies. 

The ultimate goal of ICF research is the development of an Inertial Fusion Energy Power Plant. 2 
Such a plant would most likely utilize a 2 ~ 5 MJ driver with a 500 MJ yield/target (gain >100). 
Using a 5-10 Hz shot rate, this would produce 1000 MW electrical power. With reasonable 
progress, a demonstration inertial fusion power plant could be expected to be operational in 2020 ~ 
2030. 

Cryogenic targets planned for ICF experiments are in the 
form of a liquid or solid hydrogen fuel layer which coats 
the inner surface of a low-z (e.g., polystyrene) 
microballoon. The region interior to the fuel layer (the 
“bubble”) contains low pressure hydrogen gas at the vapor 
pressure of the liquid or solid fuel layer. 

A critical concern in the fabrication of these targets is 
ensuring that the hydrogenic fuel layer maintains smooth, 
spherically symmetric uniformity with peak-to-valley 
variations on the 1000 A scale. Non-uniform or rough 
surfaces in the target or fuel can instigate Rayleigh-Taylor 
instabilities. 3 These instabilities can grow— allowing the 
cold condensed fuel to mix with the inner hot gaseous 

region and damp out the fusion reaction. To prevent Rayleigh-Taylor growth, the smoothness of 
the (plastic) shell surface must be less than 500A. 

Solid fuel layers have the advantage of being able to support the condensed hydrogen fuel layer 
naturally. However, surface smoothness is a problem. To date, no solid layering technique has 
been able to produce a smooth inner layer at the 1000 A level. Because of surface tension, liquid 
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Fig. 1. Liquid layer in a spherical shell, a 
is the radius of the vapor bubble, a x the 
amplitude of the l = 1 harmonic, b is the 
inner radius of the shell and c is the outer 
radius of the shell. 
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fuel layers will have a smooth inner surface. However, an unacceptable departure from uniformity 
arises by gravity-induced sagging of the liquid. Hence, the layer thickness will vary from thin at 
top to thick at bottom (Fig. 1). 

There are almost a dozen suggested layering methods for both liquid and solid fuel layers. For 
example fast refreeze, beta-decay heating, 4 thermally induced gradients using helium gas jets, 5 
graded intensity laser beams, 6 and microwave plasma heating. 7 The beta-layering technique is a 
volumetric heating process which relies on tritium beta-decay to produce a uniform layer, hence it 
is only good for DT layers. Thermal gradient methods induce dynamic circulating flow patterns in 
liquid fuels. The fuel can be distributed around inside the capsule by adjusting the temperature 
profile to produce a uniform thick layer. Only beta-layering has demonstrated 100 pm thick fuel 
layers. The best beta-layered surface to date has a surface roughness on 2 pm, an order of 
magnitude larger than desired. 

Plastic foams 8 can be used to trap liquid fuel. The foam density must be low (< 30 mg/cc) to 
minimize the amount of higher z materials (the carbon in the plastic). The foam cell size must also 
be small (< 1 pm) to minimize surface irregularities. If the volume of the liquid significantly 
exceeds the volume of the foam pores (overfilled), thick overfill layers will sag due to gravity. 
Thin overfilled layers (< 1 pm) may conform to irregularities commensurate with the foam cell 
size. 9 A rough estimate of the minimum overfill thickness to avoid surface imprinting (and 
Rayleigh-Taylor instabilities) would be about three times the foam cell size. ' 

Microgravity Methods 

There are two microgravity layering methods that have the potential to produce adequate liquid 
targets. These methods work by reducing the gravitational force, preventing the liquid layer from 
slumping to the bottom of the target. 

One method for reducing the gravitational force is to use the simple ffeefall method. In the absence 
of other forces, surface tension will drive the inner surface of the liquid to a spherical shape. AS 
will be shown, short range forces will center the bubble after an appropriate damping time. 
Unfortunately, for an ordinary liquid layer, the damping times are much longer than a typical ~1 
second drop time. However, if we consider an overfilled foam target with a thin ~1 pm liquid 
layer, by virtue of its smaller inertia, the thin liquid layer will approach a stable, more symmetric, 
equilibrium in a relatively short time. 

The other technique 10 ’" uses an applied magnetic field to support the liquid fuel layer against 
gravity. The method proposed here makes use of the diamagnetic nature of hydrogen. The 
magnetic dipole force per unit volume is x B* VB/Po- Balancing this against the gravitational force 
per unit volume F g = pg yields the required product of field strength with field gradient: 
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Taking x = -2.2 x 10' 6 (SI), p = 200 kg/m 3 , g = 9,8 m/sec 2 , one obtains 1.12 x 1 0 3 TVm or 1 1.2 
TVcm for deuterium, a value achievable with commercially available superconducting magnets. 

In the following we investigate the degree of symmetry attained by each of the two layering 
methods described above, i.e., the concentricity of the inner surface of the liquid fuel layer and 
examine the time-dependent motion of the liquid fuel layer. A more complete version of the 
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analysis presented here is given in ref. 1 1. We first note that the amplitudes of the perturbations 
on the inner fuel layer are extremely small compared to the radius of the target so that a 
perturbative approach is warranted. 

Assuming axial symmetry, we can then describe the radial coordinate of the inner surface of the 
liquid fuel layer, r s (0,t), by expanding it in Legendre polynomials: 

r s (0,f) =a 0 (f) + X a e (t) P<cos 0 (2) 

t- 1 

where t is time, 0 is the polar angle in a spherical polar coordinate system (Fig. 1), and is the 
perturbation amplitude of the £ tYi spherical harmonic (a e «a 0 < b < c). In the case of overfilled 
foam targets, b is replaced by d, the inner radius of the foam layer. This method of decomposition 
of surface nonuniformities (as well as absorbed laser energy and pressure nonuniformities) into 
spherical harmonics, is traditionally used in target performance and fusion gain calculations. 

We study the liquid fuel layer beginning with the energy equation: 

+ Pgm + Py + Plv) + Qvis ** 0 (3) 

with E being the kinetic energy of the whole layer, Q™ is the viscous dissipation rate, and the 
three potential energies considered here being P gm , the potential energy in the combined 
gravitational and applied B-fields or net gravitational field in the case of freefall, P y the free energy 
stored in the liquid/vapor interfacial surface, and P LV is the potential energy stored in the long 
range attractive forces between the molecules of the liquid fuel layer and the outer shell or foam 
matrix. Since all the molecules are neutral and non-polar, only the potential associated with the 
London- Van der Waals force is relevant. Each of these components can be analyzed separately. By 
using the orthogonal properties of the Legendre polynomial, all the terms in Eq. 3 can be reduced 
to functions of the mode amplitudes. 

The surface tension energy at the liquid/vapor interface (with constant surface tension y) can be 
written as: 

P y = 2nyi{e-l){t + 2){2£+l)~ l a* ' (4) 

It is noteworthy that the first £ - 1 mode results in no change of the surface area or its potential 
energy, P y . This simply means that the interior “vapor bubble” can undergo a vertical shift without 
changing its liquid/vapor surface area; accordingly, surface tension forces alone cannot 
counterbalance the gravitational forces. A similar conclusion was reached in Ref. 12 using the 
extended Young-Laplace equation of capillarity. 

The other energy terms depend on the velocity perturbations within the fluid layer. Assuming 
incompressible and irrotational flow, the velocity field is represented by v = - V § . The volume 
integrals can then be converted to surface integrals, viz: 


E = — { 
£- 2 f 

jJ <f> V<J> • n dS 

Pgm ~ ^ ^ gm V* fl dS 

(5a, b) 

ftv-j 
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v* a ds 

Qvis - T| J V v 2 - n dS 

(6a, b) 
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where n is the surface normal. 


n ■■ 
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( 7 ) 
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where P £ (|i) = dPe / d\i , is the viscosity of the liquid fuel, and O gm = U g + is the combined 
gravitational and magnetic dipole potentials per unit volume; U g =pgz, U m = \%\ B * B /2p 0 . 

The attractive interaction energy between a liquid hydrogen molecule and a shell (or foam) 
molecule is due to the long range dipole-induced dipole interaction for non-polar molecules such as 
glass or plastic. In the electrostatic (c — » °° limit), the interaction potential obeys the London-Van 
der Waals inverse sixth power law formula. However, since the dipole moments of the individual 
atoms are fluctuating at a frequency characteristic of the electronic motions (0 o , there will be a time 
between the fluctuating dipole that radiated the electromagnetic field and the reradiated field when it 
arrives back at the fluctuation and interacts with it. When the distance between two molecules, d, 
is such that the propagation time die exceeds the fluctuation period ~2 7t/co 0 , the “retardation effect” 
is significant and reduces the range of the Van der Waals forces. 13 The retardation effect can be 
quantified by multiplying the Van der Waals interaction energy between the molecules with the 
empirical expression f(p) ~ 2.45 p~ l - 2.17 p- 2 + 0.59 /r 3 , where p = o ) 0 d/c. 

t 

The total interaction energy between a single liquid molecule of mass m at position r and the 
entire plastic sphere (or foam) layer is given by the pairwise approximation. 13 With d = | r' - r \ 
this is: 


d>Lv(r) 


N eff B eff 


m 


Jv I 


volume of shell (foam) 


<*W(p) 

\r'-r | 6 


( 8 ) 


and B eff is the effective London constant defined through the weighted product for foam-filled fuel 
layers as N eff B eff = B S V S (1 -/ void ) +BHlV H / void , where N s is the number density of the solid (s) 
portion (the plastic shell) and N u is the number density of the liquid isotope (H) filling the void, 
with / void being the volume void fraction, typically f yoid = 0.95. 

The equilibrium amplitude of the first mode, a\ q , results from a balance between the net vertical 
force and the strength of the London-Van der Waals force. The layer concentricity a eq la is then 
given by 

1ST - g M(*) + E »(*) (9) 

where $'(^) is the derivative of §(x) = <^lv x 3 m c 3 1 An N eff B eff with respect to x=rlb. This 
result can also be obtained by setting the time derivatives d\ and a\ to zero in Eq. (3). 

For freefall (assuming no turbulence or vibration), surface tension drives the higher order 
equilibrium amplitudes to zero. In the case of magnetic field assisted microgravity for which B- VB 
has small non-vertical components, these higher order amplitudes, a* la, are proportional to 
a 2 y- 1 (a/L) e ' 1 where L is the scale length of the magnetic field gradient. Fortunately, the higher 
order modes are less than 0.01 % if L >5 cm and a~ c ~500|im. 
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Fig. 2a flf for hollow plastic shells: the solid line is for a 1 000 pm dia., the dashed line a 250 pm dia. 
Fig. 2b a\ q for a 100 pm thick foam shell having a density of 50 mg/cc (/ void - 0.95). 

As can be seen (Fig. 2a), a" q is directly proportional to the reduced gravitational field, eg. Its 
dependence on the thickness of the liquid layer, a‘ q ~ A 5 comes from taking the limit as A/a« 1 in 
Eq. (9). 


Fig 2b. shows that a 1 mm diameter overfilled foam target can support a liquid layer ~0.2 pm thick 
in a 1 g environment. For a given nonconcentricity, A ~ c 1/5 ; for a eq /a= 1%, A ~ 0.7 pm for e = 
1 0- 3 (free fall) and ~3 pm for e = 1 0*6 (B- VB). 
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Fig. 3 Time constants and oscillation frequencies as a function of mode number for an overfilled foam target. 


171 



Figure 3 displays the time constants as a function of liquid fuel layer thickness A for a large shell 
(c = 500 (im) with an ordinary liquid D 2 layer in conditions typical of BVB experiment (1 O' 6 g) 
and an overfilled foam target in freefall and stationary situations (10 -3 g and 1 g). The long 
damping times for the £ = 1 modes are due to the extremely weak Van der Waals 14 attractive forces 
between the liquid and the plastic shell or foam material. 

In the case of the overfilled foam layers, the decay times for all modes can be short enough to 
consider the simple freefall method, where drop times from height h, f drop = (2 h/g) l/2 , can be 
manifestly longer (nearly 1 sec for a 2 meter drop). 

Conclusions 

Symmetric cryogenic targets are critical for inertial confinement fusion. To date, no layering 
method has demonstrated the needed symmetry and surface smoothness. Microgravity techniques 
may achieve the required symmetry, but development is needed. Specifically we need to measure 
the London-Van der Waals Forces between liquid deuterium and plastic and observe the actual 
time dependent behavior of liquid layers in microgravity environments. This should lead to an 
accurate model of the liquid layer interaction that will allow the ICF community to design present 
and future ICF targets. 
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